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Abstract 

The dark energy component of the cosmic budget is represented by a self-interacting scalar 
field. The violation of the null energy condition is allowed. Hence, such component can also 
represent a phantom fluid. The model is tested using supernova type la and matter power 
spectrum data. The supernova test leads to preferred values for configurations representing the 
phantom fluid. The matter power spectrum constraints for the dark energy equation of state 
parameter are highly degenerated. In both cases, values for the equation of state parameter 
corresponding to the phantom fluid are highly admitted if no particular prior is used. 
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Observations indicate that the universe today must be in a phase of accelerated expan- 

innnn 

sion [l|, [2|, [3j, |4], |5j. In order to explain the evidences for an accelerated universe, keeping 
the traditional general relativity theory untouched, it is necessary a fluid with negative 
pressure, generally called dark energy, since it must not emit any kind of electromagnetic 
radiation. At the same time, dark energy must remain a smooth component of the cosmic 
budget, since it does not appear in the dynamics of local virialized systems, like galaxy 
and clusters of galaxies. This feature requires also a negative pressure. 

Considering a fluid with an equation of state of the type p = ap, p being the pressure 
and p the density, the parameter a must satisfy the condition a < —1/3 in order to 
drive the accelerated expansion of the universe remaining at same time an unclustered 
component of the cosmic budget: dark energy must violate the strong energy condition. 

There are some claims that the observational data favors negative values for a such 



that the null energy condition p + p > (a < — 1) is also be violated [6|. If this is the 
case, the universe may evolve towards a future singularity, since the violation of the null 
energy condition may lead to the divergence of the curvature invariants in a future finite 
proper time. This is due to the fact that the violation of the null energy condition implies 
that the density of the fluid grows as the universe expands - it is remarkable that the 
divergence in the density, and consequently in the curvature, happens in finite proper 
time and not asymptoticaly. 

It is essential to verify the strength of the evidences for a violation of the dominant 
energy condition. A lot of work has been devoted to this question, see for example [7] 
and references there in. The most simple representation of dark energy is through a fluid 
with a fixed equation of state of the type p = ap, with a < —1/3, the phantom case 
corresponding to a < — 1. In what concerns the SN la test, that is all it is needed as 
far as the equation of state of the dark energy component does not evolve. But, for the 
other tests, like those requiring perturbative analysis, the real nature of phantom field 
is crucial, that is, it is essential to know if the results are obtained considering the dark 
energy component (and in special the phantom fluid) as a self-interacting scalar field or a 
scalar field with non-canonical kinetic term^etc. There are many works exploring different 
possibilities, see for example references js, 9, 10], to quote just the more recent ones. In 
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some cases, a phantom phase can be represented in such a way that it can evolve towards 



a non-phantom configuration in the future JjJ, ll2| , leading to the avoidance of the future 



singularity. This possibility requires, of course, a non-constant a. 

In what follows we will consider this simplest case where dark energy, at least at 
background level, must satisfy an equation of state of the type p = ap, a being a constant. 
Hence, no evolution of the equation of state parameter a with time will be allowed. 
Evidently, in such situation, the precise value of a is very important in order to extrapolate 
the evolution of the universe for the very distant future. For example, if a < —5/3, it can 
appear instabilities in the perturbations in the very large wavelength limit, inducin g an 
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avoidance of big rip due to violation of the homogeneity and isotropic conditions 
Moreover, a = — 1 represents the cosmological constant, and its crucial to understand to 
which extend it is a special value in the, otherwise, continous parameter. In considering a 
constant equation of state for the dark energy fluid, the question of the future singularity 
is not addressed. 

In order to do clarify the situation concerning the evidences for a phantom cosmology, 
at least in the framework considered here, we will concentrate in two observational tests: 
the matter power spectrum and the supernova type la data. For the first case, we will 
consider the data from the 2dFGRS observational program [15J; for the second, we will 
use the data of the gold sample [4j. The model contains two fluids: one pressureless 
component, which include baryons and dark matter, and a dark energy component, rep- 
resented by a fluid which obeys the equation of state p = ap, with a in principle negative. 
The novelty of this analysis consists in the representation it will be employed for the dark 
energy fluid, a self-interacting scalar field with a specific potential, and the absence of any 
prior in the statistical analysis besides the spatial flatness. Concerning the influence of 
priors in the evaluation of cosmological parameters using observational data, in a specific 
context, see references jit], \v\. 

With the two components described above, the Einstein's equations for a flat, isotropic 
and homogenous universe described by the Friedmann-Robertson- Walker metric, reduce 
to the following expression: 

H 2 = Q m0 a- 3 + Q x0 a-^ 1+a \ (1) 



where Q m o = Qdmo + is the ratio of matter component to the critical density today, 
including dark matter (Qdmo) and baryons (fi&o)> while fl x o is dark energy ratio to the 
critical density. 

The SN la observational test constrains only the background relation, through the lumi- 
nosity distance function. However, the matter power spectrum analysis depends strongly 
on the nature of the components. For example, a fluid or a scalar field representation 
leads to complete different results for the parameter estimations. The situation is more 
delicate when components with negative pressure are considered, as it is the case for dark 
energy: a fluid description leads to an imaginary sound velocity, being unstable at small 
scales, while a field description through self-interacting scalar field implies positive sound 
velocity at sub-horizon scales [181 ]. Since the observational data for the matter power 
spectrum concerns sub-horizon modes, the specific descriptions for dark energy and dark 
matter are fundamental to interpret the observational constraints. 

To cope with the instability problem described above (which to some extend excludes 
the possibility of an ordinary fluid description for the dark energy component), the dark 
energy field will be described through a self-interacting scalar field. This is the simplest 
field description of a given component in cosmology. If it is a realistic description or a 
unique one (certainly, to some extent, it is not) is outside the aims of the present work: it 
would require, to answer this question, to know the origin of the dark energy field, what 
is object of speculation, with no clear candidate. A comparison of specific model with the 
observational data aids, of course, to shed some light on this question. 

In the absence of matter field, it is quite easy to reproduce the behaviour of a dark 
energy field through a self-interacting scalar field. The Friedmann equation coupled to a 
self-interacting scalar field and the Klein-Gordon equation for the scalar field reads, 



a' * 2 ^ 



3R =^ + V{^)a\ (2) 

(f) +3-0 = -e a , (3) 

a dcp 

where the primes mean derivative with respect to the conformal time rj defined by the 
expression dt = a(j])dr], and e = +1 is required to describe a "normal" dark energy 
fluid, while for a phantom fluid e = — 1: for a phantom field, the kinetic term must have 



the "wrong" sign. For a general equation of state p = ap, the scale factor behaves as 
a = aor] 2 ^ 1+3a \ This behaviour can be reproduced by a self-interacting scalar field with 



the form 
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V{4) = V <*V*M+ , 4>( v ) = ±2^^lln V , (4) 

1 + 6a 

Vq being a constant. It is not a surprise the appearence of an exponential potential, see 
for example reference [3]. An interesting analysis for a power law potential, implying a 
non-constant equation of state, has been performed in reference [20]. 

When matter is present, the potential (j3J) does not represent anymore exactly the dy- 
namics of the dark energy fluid. In fact, this representation is exact only in the asymptotic 
limit. When pressureless matter is present the potential that reproduces the coupled sys- 
tem dark energy/matter is more complicated, and it can not be represented, apparently, 
in a closed form using elementary functions. However, the scalar field and the potential 
for this case can be implicitly expressed in terms of the scale factor. The overall dynamics 
is accounted by the following expressions: 

r\ 2 

= a m0 a- 3 + n x0 a-^ 1+a \ (5) 

n x = ^ a' 3(1+Q) = e Y + Va\ (6) 
0' = V^V^TT^a- 3 ^)/ 2 , V = ^(1 - a)a-^ 1+a l (7) 

The relation between the scalar field and its potential to the scale factor is obtained by 
imposing that 

6' 2 

e-y + Va 2 = 8nGp x = 87iGp x0 a- (1+3a) , (8) 
6' 2 

ey-Va 2 = 8nGp x = 8nGap x0 a- {1+3a) . (9) 

When Q mQ = these expressions can be re-inserted in the Einstein's equation in order 
to have the explicit dependence of the potential in terms of the scalar field <fi. The 
cosmological constant case is reproduced in the sense that the kinetic term becomes zero 
for a = — 1, and we remain only with a constant potential term. 
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The Supernova type la analysis can be performed by using the moduli distance quantity 
defined by 

^ = 51og 10 (A/Mpc) + 25, (10) 

where the luminosity distance Dl is given by 

_c_ r dz^ 

H~oJo v/fi m0 (l + z'f + (1 - n m0 ) (1 + ^') 3(1+Q) ' 



D, 



:i + z) 



where z is the redshift. This expression is valid for a flat universe for which Q m o + Q x o = 1. 
We parametrize the Hubble parameter today writing Hq = 100 h km/ (Mpc.s). The model 
contains three free parameters, a, Qdmo and h, while the baryonic component is fixed such 
that Q b0 = 0.04. We use the SN la gold sample. The \ 2 statistics is defined by 



XSN 



N 

£ 

i=l 



o\2 



(12) 



where \x° is the observational data for the moduli distance for the i th supernova, fi\ the 
corresponding theoretical prediction and erf is the observational bar error including the 
dispersion velocity. The probability distribution function (PDF) is obtained through the 
expression 

P(h,n m0 ,a)=Ae-*s»/ 2 , (13) 

where A is a normalization constant. The PDF is three dimensional. Two dimensional 
and one dimensional PDF can be obtained integrating in one or two variables. 

Let us now turn to the perturbative analysis. We will use the Bardeen's gauge-invariant 
formalism. For the case including pressureless matter and a self-interacting scalar field, 
the perturbed equations read (see 2l|): 



V 2 $ - 3W- 



k'2 



m 1 - e— 



$" + 3 



2H' + H 2 + e- 
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k'2 



b" + 2m<p' - V 2 5(f) + e 



$ = 4nGa 2 5p + e—c 
2 



$ = ATtGa 2 5p + e-^S4 
j(f) = 40'$' - 2el> 2$ - 



+ 



(14) 

(15) 
(16) 



In these expressions, TC = a'/ a, and the subscript indicate derivative with respect to <p. 
Since the fluid represents matter, 5p = 0. The anisotropic pressure is made equal to zero 
in these equations. 



It is convenient to use the scale factor a as the new variable. The perturbed equations 
and the background relations can be re-expressed in terms of this new variable. For the 
perturbed equations we find: 
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(17) 



where A = 
definitions: 
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and the dots mean now derivative with respect to a. We have the following 



^[fi m0 + (1 - 3a)a x0 a- 3a ] 
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(19) 
(20) 
(21) 
(22) 

(23) 

(24) 



where the subscript <ft means derivative with respect to the scalar field. Moreover, k is 
the wavenumber of the perturbation coming from the Fourier decomposition and Iq = 
3.000 • hMpc is the Hubble radius today. 

We will perform a numerical integration of equations fll7lll8p . The initial conditions 



are used employing 
primordial spectrum 
in reference 



;he BBKS transfer function, and supposing a Harrison- Zeldovich 
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24j . The implementation of the initial conditions is described 



251 ]. We will compute the matter power spectrum, defined as, 
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(25) 



Sk being the Fourier component of the density contrast. As in the preceding SN case, we 
can evaluate the x 2 parameter that gives the quality of the fitting of the observational 
data by the theoretical model: 
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(26) 



where ki corresponds to the i th Fourier mode, Pj£ is the theoretical prediction for this 
mode, P^ is the corresponding observational data, and a\ its observational uncertainty. 
Since we use modes corresponding to the linear regime (scales larger than 10 Mpc), it is 
not necessary to use the full correlation matrix. 

From the Xps> we can define the probability density function (PDF) as 

P(n dm0 ,a) = Ae-*rs/2, (27) 

where A is a normalization factor. It depends, as indicated, on two free parameters, the 
dark matter fractional density Qdmo and on the equation of state parameter, a. The 
baryonic density is fixed as before. 

In figured] we display the two and one dimensional PDF when only the SN la data (gold 
sample) are used. As explained before, since the spatial section is supposed flat, there 
are three independent parameters: h, Qdmo an d a. Minimizing Xsni we ^dmo — 0.47, 
a = —2.40, h = 0.66, with Xsn = 1 11. The two (one) dimensional PDF is obtained 
integrating on the remaining one (two) parameters. The two-dimensional probability 
distributions show that values around a = —2 and h = 0.65 are favored. This is confirmed 
after marginalization: the peaks of probability occur at a = —2.29 and h = 0.66. Remark 
that the probability for a decreases after the maximum, but slowly. The extension of 
the range of a to very large negative values has increased a little the estimated value of 
h, in the direction of the value predicted by the CMB test which is around h = 0.72. 
These results agrees in their general lines with those of reference |7|], where the origin of 
the dark energy component is traced back to quantum effects, and the constraints are 
obtained by imposing that the resulting scalar field does not have a mass larger than the 
Planck mass: in this work, for the configuration corresponding to our model, there is a 
maximum of probability around a = —2. For Qdmo the analysis is a little more delicate. 
The one-dimensional PDF predicts a peak at Q m o = 0.49, a large value compared to the 
ACDM model, for which Qdmo ~ 0.25 0, 4J. But, remark that, from the two-dimensional 
PDF for a and Qdmo, it is clear that large negative values for a demands larger values for 
Qdmo- Hence, after normalization, extending the integration to deep negative values of a 
increases the predicted value for Qdmo- 
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FIG. 1: The two-dimensional PDF using the SN la data for different combinations of h, f^mo and a are 
shown in the top pannels. In the bottom panels, it is displayed the corresponding one-dimensional PDF. 




FIG. 2: The two-dimensional PDF using the matter power spectrum data for £ldm.o arid a is shown in 
the left pannels. In the center and right panels, it is displayed the corresponding one-dimensional PDF. 
Remark that the probability is almost constant. 
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FIG. 3: The two-dimensional PDF for O^mo & n d ^ using both the SN la and the matter power spectrum 
data is shown in the left pannels. In the center and right panels, it is displayed the corresponding 
one-dimensional PDF. 



The power spectrum analysis put new constraints. The power spectrum analysis 
constrains fldmo and a, the results being function of h. The best fit scenario implies 
a = —0.90, Qdmo — 0.95 with Xps = 0-38. The two and one dimensional PDF are shown 
in figure [21 It is clear that positive values for a are excluded, and after a ~ —1/3 a 
plateau is reached. Nothing special seems to happen for a = —5/3, the new critical 
point identified in references [kJ Q], but we remember that we have evaluated the power 
spectrum for the matter component. Hence, there is no contradiction with the results 



of reference 
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141 ]. After marginalization, the probability is essentially constant from 



a ~ —0.3 on; for fldmo it occurs near 1. The variation in the PDF is very small, for both 
parameter, as far as a < —0.3. 

Composing the joint PDF for both set of data we obtain the two and one-dimensional 
PDF displayed in figure [31 The maximum PDF for Qdmo is again at 0.49, the same 
position as in the SN la case. The maximum PDF for the equation of state parameter is 
at a = —2.29 as in the pure SN la case. 

The analysis here is restricted to the case where the dark energy component is described 
by a self-interacting scalar field leading to a constant equation of state. In this model, the 
evolution with a constant equation of state corresponds to a critical point in the phase 
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space, but it is not the only possibility. Introducing perturbations, the effective equation 
of state changes, and that is why instabilities do not appear even when a is negative. 
This is convenient in order to perform the power spectrum analysis. The main message 
encoded in the results obtained here, concerning the self-interacting scalar field model for 
dark energy analysed in this work, seems to be the following: there are strong evidences 
for a phantom fluid with a very negative value for the equation of state parameter a, 
mainly due to the SN la constraint; otherwise, using only matter power spectrum, the 
only clear restriction is that a must be smaller than a. ~ — 1/3. It is important to stress 
that no special prior has been used, in opposition with the analysis made, for example, 
f we particularize the value of the dark energy component for that 
3|, |4j we find essentially their results, with a peak in the probability 
distribution for a around —1. 
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